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ABSTRACT 

We calculate the time evolution of the flux, apparent size, and image centroid motion of gamma-ray burst 
(GRB) radio jets, and show that they can be resolved by the VLBA at distances of hundreds of Mpc. We find 
that GRB 030329 which showed spectroscopic evidence for an associated Type Ic supernova (SN) at a distance 
of w 800 Mpc might just be resolvable by VLBA after several months. The prospects are much better for jets 
that are oriented sideways in similar SNe with no GRB counterpart; in particular, the motion of the flux centroid 
in such jets can be detected by VLBA up to z ~ 1 even when the jet cannot be resolved. If most GRBs are 
accompanied by a SN Ib/c, then there should be a few SN/GRB jets per year within a distance < 200 Mpc and 
most of them would be oriented sideways with no 7-ray or X-ray precursor. Detection of these jets can be used 
to calibrate the fraction of all core collapse SNe that produce relativistic outflows and determine the local GRB 
rate. Overall, the rate of SNe Ib/c which do not produce a GRB at all, but rather make relativistic radio jets 
with an initial Lorentz factor of a few, may be larger by up to two orders of magnitude than the rate of those 
that produce GRBs. 

Subject headings: gamma rays: bursts — supernovae: general 



1. introduction 

Recent evidence indicates that long-duration 7-ray bursts 
(GRBs) are associated with Type Ic supernovae (SNe); of 
particular significance is the 3500 -8500 A spectrum of SN 
2003dh associated with GRB 030329 (Stanek et al. 2003), 
which was very similar to that of SN 1998bw/GRB 020405 
(Nakamura et al. 2001). This supports previous, more cir- 
cumstantial evidence, such as late time bumps in afterglow 
lightcurves (Bloom 2003) and the association of GRBs with 
central star-forming regions of galaxies (Bloom, Kulkarni, & 
Djorgovski 2002). This evidence raises two basic questions: 
(1) which fraction of all core collapse SNe produce relativistic 
outflows?, and (2) what is the probability distribution of the 
collimation angle, initial Lorentz factor, and energy output of 
these outflows? 

For every nearby (< 1 Gpc) event like GRB 030329, there 
should be hundreds of GRB jets that are not pointed at us, 
based on existing estimates for the jet opening angle (Frail et 
al. 2001). As sufficiently close GRBs can be resolved by ra- 
dio telescopes (Woods & Loeb 1999; Cen 1999; Granot, Piran 
& Sari 1999; Ayal & Piran 2001; Paczynski 2001), an effec- 
tive method to address the above two questions is to search 
for relativistic outflows in nearby core collapse SNe. In ^2] 
we calculate the expected flux, apparent size, and image cen- 
troid motion of semi-relativistic GRB radio jets viewed side- 
ways. Such jets would have no observable GRB precursor but 
can be identified and timed via their associated SN emission. 
The early images of relativistic SN/GRB radio jets may re- 
semble relativistic radio jets in quasars (Begelman, Blandford 
& Rees 1984) or micro-quasars (Mirabel & Rodriguez 1999). 
However, while quasars often inject energy continuously into 
the jet, GRB sources are impulsive. Although quasar jets re- 
main highly collimated throughout their lifetimes, GRB jets 
decelerate and expand significantly once they become non- 
relativistic, ~ 1 yr after the explosion. The hydrodynamic 
remnant of a GRB eventually becomes nearly spherical after 
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- 5 x 10 3 yr (Ayal & Piran 2001). 

2. FLUX, SIZE, AND CENTROID SHIFT, OF GRB JETS 

First, we calculate the radio flux from GRB jets observed 
from different viewing angles, 9 \, s , with respect to to the jet 
axis. We assume a double-sided jet, and calculate the emis- 
sion from both the forward jet that points towards the ob- 
server (# bs.i < 7r/2), and the opposite counter-jet ((9 Q bs,2 = 
7r _ ^obs,i > tt/2). Off-axis lightcurves from GRB jets were al- 
ready calculated using different models with various degrees 
of complexity (Granot et al. 2002), varying from simple mod- 
els (Dalai, Griest & Pruet 2002; Rossi, Lazzati & Rees 2002) 
to numerical simulations (Granot et al. 2001). Compared to 
simpler models, simulations typically show differences of or- 
der unity in the flux around or after the time of the peak in the 
light curve, a much more moderate rise in the flux before the 
peak, and a much smoother peak at # bs < 36q, where 9q is the 
initial jet opening angle. The more moderate rise before the 
peak is due to relatively slow material at the sides of the jet 
which emits much more isotropically compared to the front of 
the jet (9 < 9q) where the emission is strongly beamed away 
from off-axis observers at early times. A simple model for 
the emission from the material behind the bow shock of GRB 
jets, which essentially addresses the same emission compo- 
nent, was also investigated by Wang & Loeb (2001). 

For simplicity, we adopt an extended version of model 1 
from Granot et al. (2002); where appropriate, we mention the 
qualitative differences that are expected in more realistic jet 
models. We assume a point source that moves along the jet 
axis. Its radial location R, the lab frame time t, and the ob- 
served time f b s are related by f b S = t-Rcos9 \, s . For an on- 
axis observer (f? b s = 0) we assume a broken power law spec- 
trum (Sari, Piran & Narayan 1998). The values of the peak 
flux and break frequencies before the jet break time tj (? bsj) 
are taken from Granot & Sari (2002). 1 1> < t < ?nr, where 
'nr (fobs.NR) is the non-relativistic transition time, the tempo- 
ral scaling of the peak flux and break frequencies is modi- 
fied according to Rhoads (1999) and Sari, Piran & Halpern 
(1999). At t > ?nr the scalings are changed to those for 
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FIG. 1 . — Radio light curves (V Q bs = 43 GHz) for different viewing angles, 
and E = 10 51 erg, n = 1 cnr 3 , 8 = 0.1, e e = 0.1, e B = 0.01, p = 2.5. The 
jet that points towards us (0 o bs,l 5; 7r /2) and the counter jet that point away 
from us (0 o bs.2 = i" - Sobs.l > t/2) are shown using the same line style (the 
emission from the counter jet peaks at a later time). The SN radio emission 
is expected to be much weaker in the plotted range of observed times. 



the Sedov-Taylor regime (e.g. Frail, Waxman, & Kulkarni 
2000). The light curve for off-axis observers is then calculated 
using the appropriate transformation of the radiation field, 
^v[6Us,fob s (0ob s )] = fl 3 /v/ fl [0,f obs (0)], where a = (1 -/?)/( 1 - 
^cos^obs). In Granot et al. (2002), it was also assumed that 
fobs(0)/f o bs(#obs) = a, which is an approximation since actually 
if fobs(0)/c/fobs(#obs) = a ' an d is not very accurate for tj <t < 
fjsiR when 7 drops exponentially with radius. In this case, we 
must use the more accurate relation, f b s = t-Rcos9 b s , where 
R = J' (3(t)dt. The Lorentz factor is approximately given by 
7 « e^(t/tj)^l 2 at t < tj, 7 re fV exp(l -t/tj) at tj<t< f NR 
and 7 ss [l - (3^ R (t /t NR )~ 6 ^ 5 ] at t > f NR , where we use 

Pnr = /?(?nr) = 0.5 (7nr = 2/V3 pa 1.15) for the transition 
to the non-relativistic regime. We find 



H-(t/tj) 3 6t/l6]ct 



t < ti 



R*<{1- [exp(2f/f i -2)- l/2](6ftj/&t)}ct tj<t< ? NR , 
Rnr + [(f Anr) 2/5 - l](5/2)/3 NR Cf NR t > f NR 

(1) 

where R^ R = ^(?nr)- Following Sari et al. (1999), we use 
fobs(#obs = 0) = R/4j 2 c instead of R/I6j 2 c, and adopt their 
expression for f o bs,j(0) and the relation tj = 4f b s ,j/#o- Using 
a simple energy equation, E^o = 2E/0q = (4ir/3)nm p c 2 R 3 ~f 2 , 
and ?nrA/ = 1 — In* where x = 7nr#o, we obtain 2 
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where n = no cm" 3 is the external density and E = 10 51 £5i erg 
is the energy in the two jets. 

Figure [l] shows the radio light curves at different 9 b s , for 
both the forward and counter jets. While the optical is typ- 
ically above the peak frequency v m at tj, the radio is typi- 
cally below v m , and therefore the temporal decay slope a = 

2 f ~ 3 for 9q = 0. 1 and has only a weak logarithmic dependence on 9q. 
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FIG. 2. — The angular extent of the jet, d naK , in units of $nr (e.g. Eq.l4l. 
for the same viewing angles and jet parameters as in Fig. 



-ofln/v/ff ln? bs for 9 bs < #o changes from -1/2 to 1/3 at 
fobs,/' After f bs.m when v m sweeps past the observed frequency 
I'ob.s, a = p, while at t > t^ R a = (15p-21)/10, where p ~ 
2 — 2.5 is the power-law index of the electron energy distribu- 
tion. For 6* bs > 9o the flux still rises at f b s .j < fobs < fobs,™; for 
#obs ;$ 7r/4 the flux peaks at t \, s , p = t b s . m , while for 8 bs > tt/4, 
t P = ?nr- If fobs.p = tobs.m, the spectral slope should change from 
F u oc to i/ {x ~ p) l 2 at fobs./j, which should be easy to ob- 
serve. If t p = f NR , then fobs,,, ~ f bsj(6'obs/6'o) 2 ~ f N R^obs/ 4 for 
do < &obs <C 1, while for the counter jet f b s , p ~ ?nr, so that the 
ratio of the two peak times is ~ 4/6* 2 bs , and may be used to es- 
timate 6> obs . For 6* bs > 7r/4, t b s , p for the two jets is less than a 
factor of ~ 10 apart. The light curves for the counter jets with 



#obs,2 ^ 37r/4 (0 o bs.i < 7r/4) are all very similar, and peak at 
~ 2/?nr/c ~ 2?nr. For o b s > 9o the peak flux is a factor of a 
few larger than the flux for 6> b s = at the same f b s . For more 
realistic jet models we expect a much smoother peak and a 
somewhat smaller peak flux. We note that F^nr) oc E (Frail, 
Waxman & Kulkarni 2000), so that a larger energy implies a 
larger flux at ?nr. For example, keeping the same energy per 
solid angle, and increasing f?o from 0.1 to 0.3 (0.5), would 
increase E and therefore /v(?nr) by a factor of 9 (25), and 
would modify the afterglow light curves (Granot et al. 2002). 

Figure [2] shows the maximal angular size 3 # max , of the jets 
and counter jets from Fig. \l\ in units of 
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measured from the center. The VLBA has an angular resolu- 
tion of 4 ^170 pas at 43 GHz, and may resolve the jet around 
fNR (typically a few months to years after the SN), up to dis- 
tances of D ~ 1 Gpc. The expected peak flux at 43 GHz for a 
jet with 6* bs ~ 7r/4 at D ~ 1 Gpc is ~ 1 mJy. 

The apparent velocity of the source on the sky is v ap = 
{dOm^/ dt b & )D & Rsm0 o bJ{t - Rcos9 ob s)- For t < tj, v ap ps 
const (see inset of Fig. |3) and is ps c for 6>„b s = tt/2, sub- 
luminal for f? b s > tt/2 and super-luminal for 6> b s < 7r/2. 
For 8 bs < 1 we obtain v ap ~ 2c9 o bs/(_0 2 b& + I/87 2 ) which for 
o bs > I/7 is ~ 2c/9 bs- For 8 Q bs < I/7 we do not obtain 
the familiar result v ap ~ 7c for the afterglow image, since 

3 If both the forward and counter jets are visible, the total angular size 
would be the sum their two raax . 

4 see http://www.aoc.nrao.edu/vlba/obstatus/obssum/node30.html 
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FIG. 3. — The angle of the flux centroid relative to the center of the SN 
explosion, Scent, along the projection of the forward jet (with c9 obs < 7r/2) 
on the sky, in units of 6>nr (e.g. Eq. |4j. The inset shows decent /^obs for 
D = 100 Mpc. 



we consider a point source at a fixed angle 9 from our line 
of sight, while the edge of an afterglow image viewed on- 
axis is at 9 ~ 1 /j, where 7 changes with time (substituting 
#obs ~ I/7 in our formula reproduces this result). However, 
we are mainly interested in # ob s > #0, for which our simple 
estimate of the source size is reasonable. When 7 becomes 
;$ v a p begins to decrease. At this stage, if the jet ex- 

pands sideways significantly (i.e. 9j et ^1/7) and if the emis- 
sion from the whole jet is taken into account, then v ap ~ 7c 
and # max (?NR) ~ 9nr for all 9 b s , and not just for f9 b s > tt/4. 
Just how # ma x(?NR) changes with r? b s , depends on how close 
to spherical the jet is at ?nr. Numerical simulations show that 
the jet does not expand laterally very much before ?nr and 
may approach spherical symmetry only long after ?nr (Gra- 
not et al. 2001 ; Ayal & Piran 2001). 

Figure [3] shows the angular location of the flux centroid, 
#cent- Our results are consistent with those of Sari (1999) for 
#obs < #0 and t < ?nr- The inset shows d6 cent / Wf bs = v^/D for 
ZJ=100Mpc. Due to symmetry, 6> cent = for 9 ohs = 0,7r/2. At 
t < fNR the forward jet is much brighter than the counter jet 
for # obs ^ 7r /4, and the flux centroid largely follows the for- 
ward jet. However, when the counter jet peaks at ~ ?nr, it 
becomes somewhat brighter than the forward jet, so that the 
flux centroid gets closer to the location of the counter jet (as 
may be seen from the negative values of # cen t in Fig. 0. After 
?nr the forward and counter jets have almost the same bright- 
ness and the flux centroid moves very close to the location of 
the SN, which is midway between the two jets. This implies a 
rather large change in the location of the flux centroid ~ # NR 
for 6> b s ~ tt/4, near ?nr. For a more realistic jet model, the 
peak of the counter jet light curve is expected to be flatter and 
at a somewhat lower flux level, so that it is less clear if it will 
peak above the emission from the forward jet. Therefore it is 
not obvious whether f? cent will actually obtain negative values 
near ?nr. However, the main conclusion, that a large change 
in #cent of ~ #nr for # bs ~ 7r/4, is expected near ?nr over a 
time scale Af b s < fobs, as well as the conclusion that f? cent ap- 
proaches zero at t > ?nr, are robust (unless the two jets are not 
identical or encounter a different external density). 

The best r.m.s. error on the localization of the flux centroid 
was reported as 10 /ias (Fomalont & Kopeikin 2003). Such 
an accurate localization requires a nearby bright radio quasar 
on the sky. Compared to the best available angular resolution, 



fcent can be determined with an accuracy better by a factor 
~ 10-20. Therefore, the movement of the flux centroid on 
the sky may be detected even when the jet is not resolved (i.e. 
at early times, t <C ?nr for relatively nearby sources, or near 
fNR for more distant sources potentially up to cosmological 
distances, z ~ 1, although such sources would be dim at that 
age, ~ 0.1 mJy). 

At sufficiently late times, > ?nr ~ 1 yr, when the jets be- 
come non-relativistic and begin to approach a spherical con- 
figuration, one may estimate their physical parameters, i.e. E, 
n, and the equipartition parameters for the electrons (e e ) and 
the magnetic field (eg), in a similar way as was done for GRB 
970508 by Frail, Waxman & Kulkarni (2000). This can help 
constrain the structure of GRB jets and test if they are uniform 
or vary smoothly. 

3. APPLICATION FOR GRB 030329 

GRB 030329 was detected at a very low redshift of z = 
0.1685 (Greiner et al. 2003) or an angular sidtance D sa 
590 Mpc. Despite its low energy output in 7-rays, bumps in 
its optical afterglow light curve provide evidence for later en- 
ergy injection by refreshed shocks which bring the energy of 
the afterglow shock close to its average value for all GRBs 
(Granot, Nakar & Piran 2003). Thus we expect £51 w 1. 
Since the prompt GRB was observed, f? b s < 60 and the emis- 
sion from the counter jet should peak at ~ 2?nr ~ 2/?nr/c ~ 
5 (£51 /no) 1 1/3 yr, at a flux of ~ 20 /iJy, which would be diffi- 
cult to detect. However, f b s ,NR for the forward jet could be 
somewhat earlier, around a few months, due to light travel ef- 
fects and since the jet is still mildly relativistic at ?nr- If the 
jet spreads sideways significantly during the relativistic phase 
(0jet ~ 1/7), then its angular size after a few months should be 
~ #nr ~ 270(£ , 5i/no) I/ ' 3 M as > which just might be resolved by 
VLBA. However, if the lateral spreading of the jet during the 
relativistic stage is modest, the jet might be resolvable only 
after a few years when it becomes more spherical but rather 
dim (a few tens of /iJy). The expected shift in the flux centroid 
from early times to several months may be up to an angle of 
~ sinf? b s f?NR < 9q9nr ~ 19(£ , 5i/«o) 1//3 /ias, which might just 
be detectable with the VLBA if 9 obs ps 9 . 

4. COMPARING THE RATE OF SUPERNOVAE IB/C AND GRBS 

The rate of Type Ib/c SNe in spiral galaxies is estimated 
to be ~ 0.2 per century per 1O 1() Lb(0) (Prantzos & Boissier 
2003). The luminosity density of the local universe (Glaze- 
brook et al. 2002; Blanton et al. 2002), - 10 8 L B (0) Mpc" 3 , 
implies a rate density of SN Ib/c of ~ 2 x 10 4 Gpc~ 3 yr" 1 . The 
collimation-corrected rate of GRBs is estimated to be (Frail et 
al. 2001) ~ 250 Gpc" 3 yr" 1 . Hence, only ~ 1% of all SN Ib/c 
may be associated with GRBs. 5 

However, more SNe may have relativistic outflows with low 
Lorentz factors that would not result in GRBs (which require 
an initial Lorentz factor Tq > 100), but rather in UV (for 
Tq ~ 10) or radio (for To < 3) transients only. The obser- 
vational constraints on the rates of such transients are weak. 
Calibration of the statistics of relativistic radio jets in core 
collapse SNe can be used to infer the rate of such transients 
(which should occur on the rare occasions when the same jets 
are viewed on-axis). It can also provide new and more reliable 

5 Norris (2002) has a more optimistic prediction that > 25% of all SNe 
Ib/c produce a sub-class of low-luminosity GRBs similar to GRB 980425/SN 
1998bw. 
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evidence for the collimation of GRB jets, and an independent 
estimate for the distribution of the collimation angles. 

5. CONCLUSIONS 

We have calculated the radio light curves, and the evolution 
of the apparent size and flux centroid (FC) of GRB jets viewed 
sideways. As the jets do not point at us, they will have no 7- 
ray precursor, but will instead be preceded by a Type Ib/c SN. 
A ~ 1 yr old GRB remnant at D ~ 100 Mpc, is predicted to 
have a characteristic radio flux of ~ 100 mJy, an image size of 
~ 1 mas, and FC motion of ~ 20 /xas week" 1 . Such a source 
can be resolved by VLBA at D < 1 Gpc, while the motion of 
its FC might be monitored up to z ~ 1 (although at z ~ 1 it 
would be very dim, ~ 0.1 mJy). The apparent size of the jet 
or super-luminal motion of its FC within the first few months 
after the SN, may provide evidence for relativistic motion. 

A relativistic jet of length ct may also serve as a yardstick 
for constraining cosmological parameters. However, the re- 
quired precision for this purpose may not be attainable if the 
jet orientation is not well known or the surrounding medium 
is inhomogeneous. 

For an off-axis jet, there should be a time dependent linear 
polarization, which peaks near the time of the peak in the light 
curve and slowly decreases with time as the jet becomes more 
spherical and symmetric around the line of sight (Granot et al. 
2002). If the jet is resolved, then polarization maps could be 
generated, as is commonly done for extragalactic radio jets. 
This could reveal the magnetic field geometry and orientation 
in the jet, and whether it has a large scale ordered component 
(Granot & Konigl 2003). 

The existence of an early phase during which the emis- 
sion of the jet peaked in the UV (7 > 10) can in principle 
be inferred from the ionization cones preceding the jet in the 
surrounding gas (Perna & Loeb 1998). Since the recombina- 
tion time of the gas is ~ 10 5 Hq' yr, these cones should exist 
for long times after the SN explosion. However, the sepa- 
ration between the ionization fronts and the edges of the jet 
grows large only after the jet becomes non-relativistic. At 



these late times, one may detect emission lines from highly 
ionized, metal rich gas (Perna, Raymond, & Loeb 2000) that 
reflect the hardness of the emission spectrum of the jet at ear- 
lier times when it was highly relativistic. Detection of ioniza- 
tion cones can be used to infer the early opening angle and 
spectral flux of the jet at different frequencies (corresponding 
to the the ionization state of different ions). The latter can be 
used to estimate the initial Lorentz factor of the jet, To in the 
range ~ 5-20, or determine if To < 5 (no ionization cones), 
or To > 20 (ionizing extending up to soft X-rays). 

There is strong evidence connecting GRBs with SNe Ic that 
have a large kinetic energy, > 10 52 erg (termed 'hypernovae' 
by Paczynski 1998), and have a distinct spectral signature, as 
was observed for SN 1998bw and SN 2003dh. The search 
for GRB radio jets in SN with such a spectrum is particularly 
interesting, as it could show whether all such SNe produce 
GRBs. For 6* b s < #0 the existence of relativistic jets can be 
revealed by their high brightness temperatures (^> 10 1 1 K, e.g. 
Kulkarni et al. 1998; Li & Chevalier 1999). The fraction of 
SN Ib/c that produce GRB jets can help determine the local 
GRB rate and the distribution of 9q. 

It may prove interesting to search for a correlation between 
the value of To as estimated from the ionization cones and 
the spectrum of the SN. A correlation with ejecta energy and 
abundance patterns (Maeda & Nomota 2003) may show a 
continuous change in the SN spectrum as a function of To 
of the bipolar jets, which might indicate that these jets are in- 
timately related to, or perhaps are, the main cause of the SN 
explosion. 



We thank Josh Winn, Mark Reid and Re'em Sari for useful 
discussions. This work was supported in part by the Institute 
for Advanced Study, funds for natural sciences (J.G.) and NSF 
grants AST-0071019, AST-0204514 & NASA grant NAG 5- 
13292 (A.L.). A.L. acknowledges support from the IAS at 
Princeton and the J.S. Guggenheim Memorial Fellowship. 



REFERENCES 



Ayal, S., & Piran, T. 2001, ApJ, 555, 23 

Begelman, M. C, Blandford, R. D., & Rees, M. J. 1984, Rev. Mod. Phys., 
56, 255 

Blanton, M. et al. 2002, ApJ, submitted, astro-ph/0210215 

Bloom, J. 2003, to appear in the Proceedings of "Gamma Ray Bursts in the 

Afterglow Era, Third Workshop" (Rome, Sep. 2002); astro-ph/0303478 
Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2002, AJ, 123, 1111 
Cen, R. 1999, ApJ, 524, L51 

Dalai, N., Griest, K., & Pruet, J. 2002, ApJ, 564, 209 

Fomalont, E. B., & Kopeikin, S. M. 2003, astro-ph/0302294 

Frail, D. A., Waxman, E., & Kulkarni, S. R. 2000, ApJ, 537, 191 

Frail, D. A. et al. 2001, ApJ, 562, L55 

Glazebrook, K. et al. 2003, ApJ, 587, 55 

Granot, J., Piran, T., & Sari, R. 1999, ApJ, 513, 679 

Granot, J., et al. 2001, in 'GRBs in the Afterglow Era', ed. E. Costa, F. 

Frontera, & J. Hjorth (Berlin: Springer), 312 
Granot, J., Panaitescu, A., Kumar, P., & Woosley, S. E. 2002, ApJ, 570, L61 
Granot, J., & Sari, R. 2002, ApJ, 568, 820 

Granot, J., & Konigl, A. 2003, ApJ, 594, L83 

Granot, J., Nakar, E., & Piran, T. 2003, preprint 1 astro-ph/0304563 1 
Greiner, J., et al. 2003, GCN circ. 2020 
Kulkarni, S. R., et al. 1998, Nature, 395, 663 



Li, Z. Y., & Chevalier, R. A. 1999, ApJ, 526, 716 

Maeda, K., & Nomoto, K. 2003, ApJ, submitted 1 astro-ph/0304172 1 
Mirabel, I. F. & Rodriguez, L. F. 1999, ARA&A, 37, 409 
Nakamura, T., Mazzali, P. A., Nomoto, K„ & Iwamoto, K. 2001, ApJ, 550, 
991 

Norris, J. P. 2002, ApJ, 579, 386 
Paczynski, B. 1998, ApJ, 494, L45 
Paczynski, B. 2001, Acta Astronomica, 51, 1 
Perna, R. & Loeb, A. 1998a, ApJ, 501, 467 

Perna, R., Raymond, J., & Loeb, A. 2000, ApJ, 5 33, 658 

Prantzos, N., & Boissier, S. 2003, A&A in press 1 astro-ph/0305376 1 
Rhoads, J. E. 1999, ApJ, 525, 737 

Rossi, E., Lazzafi, D., & Rees, M. J. 2002, MNRAS, 332, 945 
Sari, R. 1999, ApJ, 524, L43 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 
Sari, R., Piran, T., & Halpern, J. P. 1999, A pJ, 519, L17 
Stanek, K. Z., et al. 2003, ApJL, submitted; astro-ph/0304173 
Wang, X., & Loeb, A. 2001, ApJ, 55 2, 49 
Woods, E., & Loeb, A. 1999, astro-ph/9907110 



